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Abstract 

We study the flavour changing neutral Higgs boson decays that can be in- 
duced from genuine supersymmetric particles at the one-loop level and within 
the context of the Minimal Supersymmetric Standard Model. We consider all 
the possible flavour changing decay channels of the three neutral Higgs bosons 
into second and third generation quarks, and focus on the Supersymmetric- 
QCD corrections from squark-gluino loops which are expected to provide 
the dominant contributions. We assume here the more general hypothe- 
sis for flavour mixing, where there is misalignment between the quark and 
squark sectors, leading to a flavour non-diagonal squark mass matrix. The 
form factors involved, and the corresponding Higgs partial decay widths and 
branching ratios, are computed both analytically and numerically, and their 
behaviour with the parameters of the Minimal Supersymmetric Standard 
Model and with the squark mass mixing are analyzed in full detail. The 
large rates found, are explained in terms of the non-decoupling behaviour of 
these squark-gluino loop corrections in the scenario with very large super- 
symmetric mass parameters. Our results show that if these decays are seen 
in future colliders they could provide clear indirect signals of supersymmetry. 
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1 Introduction 



The indirect searches of supersymmetric (SUSY) particles via their radiative 
effects into low energy observables have received a lot of attention in the 
last decades The main motivation for these searches is that they can 

provide valuable clues in the way towards the final discovery of supersym- 
metry, even in the most pessimistic scenario where the SUSY spectrum is 
too heavy as to be directly produced in the present or forthcoming colliders. 
In this concern, the flavour-changing neutral current (FCNC) processes are 
'ideal laboratories' where to look for these indirect SUSY signals or any other 
radiative effects from new physics beyond the Standard Model (SM) of par- 
ticle interactions, since the SM predicts negligible rates for these processes. 
We are interested here in the FCNC effects that are generated from genuine 
SUSY radiative corrections within the context of the Minimal Supersymmet- 
ric Standard Model (MSSM) Our interest is focused particularly on the 
SUSY radiative effects that induce scalar flavour changing interactions [Q. 

The strong suppression of any FCNC process within the SM is due to the 
absence of tree-level flavour changing (FC) interactions and the GIM cancel- 
lation mechanism ||^ that operates beyond the tree-level. Within the MSSM, 
the scalar FC interactions are also absent at tree-level, but they can be gener- 
ated quite efficiently at the one-loop level and lead to sizeable contributions 
at some regions of the MSSM parameter space. This is mainly because 
the GIM suppression mechanism does not necessarily operate in the genuine 
SUSY radiative corrections . Special mention deserves the SUSY one-loop 
radiative effects that modify the tree-level relations between the down-type 
fermion mass matrices and their corresponding Yukawa coupling matrices 0, 
since they can induce FC Yukawa interactions if there is change of flavour in 
the internal squark lines This has been the subject of numerous studies 
in the last years, because these SUSY radiative effects are very significant at 
large tan (3 values . This tan /? reassures the ratio of the vacuum expec- 

tation values of the two Higgs doublets and is predicted to be large in some 
SUSY-GUT models with top-bottom Yukawa coupling unification, being of 
order mt/rrih ~ 50 [|r^. On the other hand, the phenomenological interest of 



these tan P enhanced SUSY radiative corrections leading to FC neutral Higgs 
interactions is that they affect a wide range of low energy mesonic processes 
and can be in conflict with present experimental data, including — 
mixing P,[Tl|,|l2|,|l^, leptonic B meson decays |P,|l4|,[l^,|l6|,0,|l8|], as well as 
some other CP-conserving |1^,|2^ and CP- violating observables related to 
the K and B meson physics [^]. Similar enhanced tan/5 effects have 
been found at 6 — S7 decays [^] as well. 

Another studies at higher energy processes of FC scalar interactions be- 
ing induced from SUSY loops have also been performed in the literature. 
These include rare Z boson decays [p4| , p5| and rare top decays p6| , p7| , p8| , |29 



both being significant at particular regions of the MSSM parameter space. 
Specially large are the branching ratios for the Z ^ hs decays, which can 



2 



reach 10~® for large tan/3 values and for the top rare decay into neutral 
MSSM Higgs bosons, t cH, which can be as large as 10~^ [^^, and both 
rates are largely dominated by the SUSY radiative effects from squark-gluino 
loops. With such a large rates, these Z ^ bs decays could be accessible at a 



future Giga-Z collider p5[, and the t cH decays could be observed at the 



forthcoming CERN-LHC collider [|g . 

Our subject of study here is the flavour changing neutral Higgs boson 
decays (FCHD) that can be induced from genuine SUSY particles at the 
one-loop level. This is a subject closely related to the previous ones but, to 
our knowledge, has not been analyzed in the literature yet. We will work 
here in the MSSM context and consider all the possible FC decay channels 
of the three neutral MSSM Higgs bosons into second and third generation 
quarks, namely, H — >• bs, sb, tc, ct, with H = ho, H^, and, as usual, ho 
and Ho are the lightest and heaviest CP-even bosons, respectively, and Ao is 
the CP-odd boson. Our study is devoted to the second and third generation 
quarks because the squark mixing between these two generations, which is 
been assumed here to induce these decays, is the less constrained experi- 
mentally 0. We have focused on the SUSY radiative corrections from the 
SUSY-QCD sector which, by analogy to the previous studies, are expected 
to provide the dominant contributions to the FCHD rates, as their strength 
is driven by the strong coupling constant. 

More concretely, we will compute here the one-loop contributions of 
0{oPg) to the partial decay widths of the neutral Higgs bosons coming from 
loops of gluinos and third and second generation squarks. We will present 
the analytical results from our diagrammatic computation in terms of the 
involved form factors and we will analyze in full detail the numerical values 
of the Higgs partial widths as a function of various relevant MSSM param- 
eters. Basically, the yU mass parameter, the gluino mass, the squark masses, 
the pseudoscalar mass and tan 

Regarding our hypothesis on the generation of squark flavour mixing, we 
consider here the most general scenario, which is called in the literature non- 
minimal flavour scenario. This occurs when the squark mass matrices are 
not flavour- diagonal in the same basis as the quark mass matrices. Thus, 
when the squark mass matrices are diagonalized, FC gluino-quark-squark 
couplings arise for the mass eigenstates, and these induce in turn the FCHD 
via squark-gluino loops, which are the subject of our study. It is worth to 
recall that this hypothesis of misalignment between the squark and quark 
mass matrices is present in the most general parameterization of the MSSM 
and can lead to dangerous FC effects in conflict with experiment. Specially, 
the data on i^^ — and — mixing impose severe constraints on 
the mixing involving the first generation This is why we focus on the 
mixing between the second and third squark generations, which is practically 
unconstrained 

Attempts to solve the previously commented SUSY flavour problem of 
the MSSM include flavour-blind SUSY breaking scenarios, as in Minimal 
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Supergravity, in which the sfermion mass matrices are flavour diagonal in 
the same basis as the quark mass matrices at the SUSY-breaking scale, but a 
small amount of non-minimal flavour mixing is generated as the squark and 
quark masses are evolved down, via renormalization group equations, to the 
electroweak scale. The actual estimates of these radiatively induced flavour 
off-diagonal squark squared mass terms ||31|,|3^ indicate that the largest ones 
are those referred to the SUSY partners of the left-handed quarks. All, 
since these scale with the squared of the soft SUSY breaking masses, in 
contrast to the Alh (or A^l) and Arji terms that scale with one or zero 
powers respectively of the soft SUSY breaking masses. Thus, the hierarchy 
All >> Alr >> Arr is usually assumed. These same estimates also 
indicate that the All terms generating the mixing between the second and 
third generation squarks can be numerically significant because of the third 
generation quark mass factors involved. 

In our analysis of the FCHD we will assume the simplest and theoretically 
better motivated hypothesis, where the only non-zero off-diagonal squark 
squared mass entries in the d-sector and u-sector are for and citL mix- 
ing, respectively. We will parameterize these flavour off-diagonal entries, a 
la Sher |Q, simply by A^^ = XllMqMq,, where Mq and Mq, are the two 
corresponding involved soft SUSY breaking masses. The dimensionless pa- 
rameter Xll is considered here to be the only free-parameter characterizing 
the flavour mixing strength and, for the numerical evaluations, it will be 
taken in the conservative range, < A^^ < 1, which is perfectly allowed by 
the present data (for a summary of present bounds on the Xlr, Xrl, ^ll and 
Xrr parameters, see 

The paper is organized as follows. We present in section ^ the flavour 
non-diagonal squark squared mass matrices and write the relevant quark- 
squark-gluino interaction terms that are generated after rotating to the mass 
eigenstate basis. The computation of the one-loop corrections to the form 
factors and FCHD widths and their analytical results are presented in sec- 
tion ^. The numerical analysis of the FCHD rates and a detailed discussion 
on their dependence with the MSSM parameters and with the Xll parameter 
are included in section ^. After scanning the MSSM parameter space, par- 
ticular regions where the FCHD rates are maximal are detected. The values 
of the FCHD branching ratios as a function of the A ll parameter are then 
analyzed at these regions. We devote section ^ to the study of the SUSY de- 
coupling properties in these FC observables. After performing a large SUSY 
mass expansion, we will show in this section that the SUSY radiative effects 
from squark-gluino loops indeed do not decouple in the FCHD. This explains 
why the numerical size of the FCHD rates found in this work are so large. 
Finally, the last section is devoted to the summary and conclusions. 
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2 Flavour Changing in the SUSY-QCD sec- 
tor of the MSSM 

As recalled in the introduction, in the MSSM there is a new source of flavour 
changing phenomena coming from a possible misalignment between the rota- 
tion that diagonalizes the quark and squark sectors. The squark mass matrix, 
expressed in the basis in which the squarks fields are rotated parallel to the 
quarks (called the super-CKM basis), is in general non-diagonal in flavour 
space, and provides new FC effects. 

Assuming that FC squark mixing is significant only in transitions between 
the third and second generation squarks, and that there is only LL mixing, 
given by the Sher type ansatz |3^, the squark squared mass matrices in the 
{cL,CR,iL,ifi) and {sL,SR,bL,bji) basis, respectively, can be written as follows. 
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(2.1) 



(2.2) 



(2.3) 



ruq, T^, Qq are the mass, isospin and electric charge respectively of the quark 
q; mz is the Z boson mass, and sw is the sine of the weak angle 6*14^. The 
involved MSSM parameters in the SUSY-QCD sector usual, the gluino 

mass, Mg, the /i-parameter, the soft SUSY breaking masses Mq, Mjy, Mjj 
and the soft SUSY breaking trilinear parameter, A. Notice that we have used 
in eqs.(|2.1|), (|2.2|) and (p.3|) a self-explanatory notation for the flavour indices. 



and that, due to the SU(2)l invariance, Mq ^ = Mq ,, and Mq ^ = Mq ^^. 

In our previous parameterization of flavour mixing in the squark sector, 
there is only one free-parameter, the Xn parameter, that characterizes the 
flavour mixing strength. Notice that, for the sake of simplicity, we are assum- 
ing the same notation for the Xll parameter in the t — c and h — s sectors, 
which from now on will be called simply A. For the numerical estimates 
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in the following sections, the A values will be taken in the moderate range 
< A < 1. Obviously, the choice A = represents the zero-flavour mixing 
case. 

In order to diagonalize the two previous 4x4 squark mass matrices, two 
rotation 4x4 matrices, R^^^ and R^'^\ one for the up-tjpe squarks and the 
other one for the down-type squarks, are needed. Thus, the squark mass 
eigenstates, qa, and the interaction squark eigenstates, q'^, are related by. 



(2.4) 
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(2.5) 



Once the squark mass matrices are diagonalized, one obtains the squark 
mass eigenvalues and eigenstates that obviously depend on A. Again, the 
choice A = represents the case where no flavour mixing occurs, and recovers 
the usual pairs of physical squarks, (61,62), (si,52), (^1,^2), and (£1,02) with 
the usual mass patterns. 

For illustrative purposes, we show in fig.|l| the physical masses of the down- 
type squarks, Mj^^^^, as a function of A, where, for definiteness, we have 
chosen all the SUSY mass parameters equal to 1 TeV and tan (3 = 40. We can 
see clearly the numerical differences among the four squark masses M^j 2 3 ^ 



(we take the convention Mg^ 



> Mg^ > Mqg > Mg^) and the corresponding 



ones for A = that do not include any flavour changing. In the d sector, there 
are two eigenvalues, Mj^ and Mj^ that do not change much as we increase 
the value of A. These correspond, when A = 0, to Si and §2 respectively. On 
the other hand, Mj^ increases with A while Mj^ decreases. These last ones 
correspond, when A = 0, to 61 and 62 respectively. 

The numerical analysis for the up-type squarks is completely analogous 
to the previous one of the down-type squarks and we dont show it here 
explicitly. Concerning the mass pattern and its dependence with A, there 
are two eigenvalues, M^^ and M^^ that do not change much as we increase 
the value of A, and that correspond, when A = 0, to Ci and £2 respectively. 
On the other hand, M^^ increases with A while M^^ decreases and these last 
ones correspond, when A = 0, to ti and ^2 respectively. 

Finally, notice that, for the numerical analysis of the FCHD rates in the 
next sections, only A values in the < A < 1 range that lead to physical 
squark masses above 150 GeV will be considered. The present experimen- 
tal lower mass bounds on the squark masses of the first and second squark 



generation are actually even more stringent that this value |3^], but we have 
chosen here this common value of 150 GeV for simplicity and definiteness. 
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Figure 1: Physical masses of the down- type squarks as a function of A when 
all the SUSY mass parameters are equal to 1 TeV and tan f3 = 40. The solid 
line corresponds to Mj^ , the long dashed line corresponds to M^^ and the 
other two correspond to Mj^ and M^^ (short-dashed and semi-short-dashed) 
respectively. 



Similarly, in view of the present experimental lower bounds on the gluino 
mass [0, we will consider here Mg values above 200 GeV. 

The previously introduced intergenerational mixing effects in the squark 
sector lead to strong FC effects in processes involving neutral currents through 
the quark-squark-gluino interaction Lagrangian, which can now be written 
in the squark mass eigenstate basis as, 

^9, Q, q) = -V2gst (^R^^J*gujL + Rfi*gd*J)L - R^^]*gujR - Rfi*gd*J)R 

+R^^J*gU*^CL + Rfa*9d*aSL - R^2a*9KcR " + k.C 

(2.6) 

with a = 1,2,3,4 and we have used the short notation tgq*q = t1-g^q*'^q^ , 
where t"^ are the standard S'f/(3)c generators. For simplicity, we will omit 
the color indices from now on. 

In this way, if there is misalignment between quark and squarks, new FC 
effects will appear in neutral currents processes and, in particular, in the 
neutral Higgs decays of the MSSM into quarks that we are interested in. 
These will occur, mainly, via loops of squarks and gluinos and through the 
flavour non-diagonal q—q—g couplings of eq.(p.6|) which in turn have emerged 



from the flavour non-diagonal squark squared mass matrices of eqs.(2.1) and 



(|2.2| ). These effects will be driven by the strong coupling constant a^, and 
therefore they are expected to be numerically large. We will study those 
effects in full detail in the forthcoming sections ^, ^ and |^. 
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3 Generating flavour changing Higgs decays 
from squark-gluino loops 

We present in this section the computation of the one-loop radiative cor- 
rections to the FCHD of the neutral MSSM Higgs bosons into second and 
third generation quarks, H — > bs, sb,tc,ci. Notice that, the present the- 
oretical upper mass bound nih^ < 135 GeV makes the lightest Higgs 
boson ho unable to decay into tc or ct. Therefore, just the following channels 
will be considered. Ha — > bs, sb with Ha = ho. Ho, Ao and Ha — > tc, ct with 
Ha = Ho, Ao- We will focus on the radiative corrections that come from 
the SUSY-QCD sector, and more specifically those from loops of squarks of 
type s, c, b and i, and gluinos g. These will provide contributions to the FC 
partial widths of and are expected to be the dominant ones. 

Notice that, the computation of these FC partial widths are relatively 
easy and do not require renormalization, because these decay processes do not 
proceed at the tree level in the MSSM. One just computes the different one- 
loop diagrams that contribute to the process and the final result obtained, 
after adding up all of them, must be finite, since no lowest order interaction 
could absorb the left over infinities. 

In order to present our results for the partial decay widths, it is convenient 
to define the one-loop effective interaction term associated to each decay 
H qq' in the following compact form: 



iF 



-iguM{Fl' {H)Pl + F^" {H)Pn)vAP2)H{ps), (3.1) 



where g is the SU{2)i gauge coupling and F^ and are the form factors 
associated to each chirality projection L and R respectively. The values 
of Fl and F^ follow from the explicit calculation of the one-loop vertices 
and mixed self-energies and are of order as- These one-loop diagrams from 
the SUSY-QCD sector are depicted in fig.|^ for the H ^ bs decays, with 
H = ho. Ho, Ao- The corresponding diagramas for the tc decays are identical 
but replacing d —>■ u, b ^ t and s — > c. 



d.3 



H 



9a I, 




H 




Figure 2: One-loop diagrams from the SUSY-QCD sector for the decay H — *■ 
bs 

For the mixed self-energy diagrams it is useful to define the following 
decomposition. 
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(3.2) 



and similarly for S*'^(/c) by replacing the quark flavour indices correspond- 
ingly, namely, h ^ t and s — > c. 

By computing the diagrams in flg.^, we get the following results for the 
d-sector form factors in terms of the scalar one- loop integral functions Bq, 
Bi, Co, Cu and C12, as deflned in refs. PB| , 
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Similarly, by computing the equivalent diagrams for the tc decays, we flnd 
the expresions for the w-sector form factors, Fj^{Ha) and Fj^{Ha), which are 
identical to the previous results for F^^{Ha) and F^{Ha) respectively, but 
replacing the flavour indices correspondingly, namely, b ^ t, s ^ c and 
d —>■ u. 

In the previous formulas, the gHagcgp the Higgs-squark-squark cou- 
plings in the mass eigenstate basis, which are collected in the Appendix 
A. The Higgs-quark-quark couplings are given by, Gnqq = — ^ J^"!!^ ^ , for 

s,b, and Gnqg = - ,^ZZlm q 



2Mw cos (3 ' 

2Mivsm/3' """" ^ '^^^ factors depend on 

the particular Higgs decay. These are, = (cos a, sin a, i cos /5), k'I = 



sin a, cos a, i sin /5) , and k)j 



{u4) 



K 



{u,d)* 



for the decays of Ha = {ho, Ho, Ao) 



respectively. The last parenthesis in the second line of all the form factors 
refer to the arguments of the one-loop Co, Cu and C12 functions. Finally, the 
one-loop contributions from the SUSY-QCD sector to the mixed self-energies 
appearing in the previous equations are the following, 

\k') = -'^-^B.ik^MlMl.jRtfRtf* 



Y^tc,bs 1 1 2\ 



37r 



(3.4) 
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2a. 



m,^Tik') = -^M,Bo{e,MlMlJRt'^Rt'^* (3.5) 

The expressions for the form factors involved in the Ha — > ct, tc decays have 
been checked to be in agreement with the previous resuhs obtained in ||2^ 
for the t cHa decays. 

We next present the resuhs of the partial decay widths for the decays 
Ha — > tc, ci where Ha = Ho, Ao and Ha bs, sb where Ha = ho,Ho,Ao, in 
terms of the above form factors. Since we are assuming that the final states 
qq' and q'q are not experimentally distinguishable, the final results for the 
total partial widths are got by adding the two corresponding partial widths, 
and this will be denoted from now on by r{H —>■ qq' + q'q). These results are 
as follows. 
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(3.7) 



where A 2 {l,x ,y ) = y [1 — (s + yY] [1 — (s — yY]. Notice again that, due 
to phase space, the lightest Higgs boson ho cannot decay into tc or ci. 

Some comments are in order. First, we have checked explicitly that the 
previous results for the form factors are finite, as expected. Second, we 
can see that the parameter A, which characterizes the FC, does not appear 
explicitly in eqs.(|3.3|) and (p.4|), but it appears implicitly. As we have already 
explained, the effect of FC is due to misalignment between the quark and 
squark mass matrices and is parametrized through non-diagonal terms in 
the squark squared mass matrices containing the parameter A. Therefore, 
the dependence on A is hidden in the respective rotation matrices used to 



diagonalize the squark squared mass matrices, R^^l and R^aSi obviously 



in the physical squark mass values, which appear as arguments of the scalar 
one-loop integral functions obtained in the one-loop calculation. 
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On the other hand, the MSSM parameters tan /3, /i and Mg will be crucial 
for our phenomenological analysis of section ^ where we will study numer- 
ically the size of these loop induced FCHD as a function of these MSSM 
parameters. The dependence on /i and tan j3 in the form factors, and hence, 
in the decay widths, is complicated and is hidden inside couplings and inside 
the squark mass matrices. Mg appears as an argument of the scalar one-loop 
integral functions, and also it appears explicitly as a multiplicative factor in 
both the vertex Co contributions and the scalar parts of the mixed self-energy 
diagrams. The global dependence on the angle (3 is more complicated, as it 
emerges in most of the terms. The explicit (3 dependence through LR flavour 
preserving squark mixing, namely, ^(ylfc tan /5) and mt^c{.^t,c— ^^ cot 
shows that the contribution from these mixing terms in the down sector grows 
with tan/5 and its dependence on A;, ,,, for large values of tan/5, is expected 
to be mild, while in the up sector, it is less tan/3 dependent for tan/5 > 1 
and the trilinear term, c, can acquire more importance. 

At the end, one crucial parameter for our complete analysis will be /i, 
which appears in the squark mass matrix term in the way exposed above, 
but also, directly in the Higgs-squark-squark couplings, as can be seen in 
the Appendix A. As we will show in the next sections, large values of can 
produce sizable contributions to FCHD. 

All these complicated dependences with the different MSSM parameters 
will be clarified in the analysis performed in the next section, and the an- 
alytical behaviour will be explained in more detail in the large SUSY mass 
limit performed in section |^, where we study the behaviour of these PC de- 
cay processes in the scenario where a very heavy supersymmetric spectrum 
is considered. 

4 Numerical Analysis of the FCHD rates 

In this section we estimate numerically the size of the loop induced FCHD 
as a function of the MSSM parameters and the A parameter. The MSSM 
parameters needed to fully characterize and evaluate the T{Ha hs + sb) 
and T{Ha —>■ tc + ct) partial widths, in our simplified scenario, consist of 
the following six parameters: m^, tan/5, ^, Mg, Mo, and A, where we have 
chosen, for simplicity. Mo as a common value for the soft squark breaking 
mass parameters. Mo = Mq ^ = Mj^ ^^ j^ = M^ j-^ j,-,, and all the various 
trilinear parameters to be equal At = Ab = Ac = Ag = A. These parameters 
will be varied over a broad range, subject only to our requirement that all the 
squark masses be heavier than 150 GeV, and that the gluino mass be heavier 
than 200 GeV. On the other hand, the extra parameter. A, which is the only 
one measuring the PC strength, will be varied in the range < A < 1, and 
only values leading to Mg > 150 GeV will be considered. 

In order to estimate the FCHD rates over the MSSM parameter space 
and the A interval chosen, and with the motivation in mind of finding the 
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r(h<,-+s b+bT> r(ho s b +b"s) 




Figure 3: T{ho — > + sb) in GeV as a function of (tan/5, rriA (GeV)) 
(a), {fi (GeV), tan/3) (b), (/i (Gel^), (Gel^)) (c) and (M^ (Ge\/), A 
(GeV)) (d). The regions of the parameter space not plotted are the ones 
that give non allowed values for the squark masses {Mq < 150 GeV^). The 
values of the different MSSM parameters that have to be fixed in each plot 
have been chosen correspondingly to be: /i = 1500 GeV, Mo = 600 GeV, 
Mg = 300 GeV, A = 200 GeV, tua = 250 GeV, tan/5 = 35 and A = 0.5. 

region of the MSSM parameter space where these rates are maximal, we will 
proceed as follows. We will first fix the A parameter to one concrete value 
in the range < A < 1 and vary the rest of parameters over the whole 
allowed range in order to detect such maximizing region. Next, we will fix 
all the MSSM parameters to some specific values belonging to this region 
and study the widths and branching ratios as a function of the A parameter. 
This will give us finally the maximum size of these observables as a function 
of A as well, and if any planned experiment is able to observe and measure 
these FCHD rates with good accuracy, their experimental value could serve 
to extract either the preferred A value or, in the worst case, an upper bound 
on it. 

The figs.^ through 1^ display the numerical results for the T{Ha —>■ bs + sb) 
and T{Ha —>■ tc + at) partial widths as a function of the six previous MSSM 
parameters, uia, tan/?, fj,, Mg, Mq and A, and for the specific value A = 0.5. 
In these figures, the MSSM parameters have been grouped in pairs in order 
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rCH^ -» s b + br|i r(H„^£b+bT) 




Figure 4: Same as in fig.|] but for T{Ho ^ bs + sb). 



to illustrate simultaneously the behaviour of the FCHD widths with both 
parameters. These pairs of parameters have been chosen to be: (tan/3, rriA), 
(/i, tan/5), (/i, Mg) and (Mq, A), and correspond to the plots (a) (upper left 
panel), (b) (upper right panel), (c) (lower left panel) and (d) (lower right 
panel), respectively, in these figures. We have made these four plots for each 
neutral Higgs boson. figs.|^, ^ and |^ correspond to T{ho ^ bs + sb), T{Ho 
bs + sb) and T{Ao ^ bs + sb), respectively, and figs.| and |^ correspond to 
T{Ho — > tc + ct) and T{Ao tc + ci), respectively. Notice again that the 
lightest Higgs cannot decay to tc nor ct, due to phase space. 

In the following we discuss separately the T{Ha ^ bs + sb) and T{Ha — *■ 
tc + ct) cases due to their different behaviour with some of the MSSM pa- 
rameters. 

First, we will focus on T{Ha ^ bs + sb) (see figs.|], | and H, for the ho, 
Ho and Ao respectively). Let us start with the first plot, (a), of these figures 
where we study the combined behaviour with tan/5 and m^. The first clear 
behaviour of these decay widths is the fast growing with tan (3. Notice that 
it is the same behaviour for the three neutral Higgs bosons, despite the fact 
that the growing with tan (3 is better seen in the case of the Ho and Ao than in 
the lightest Higgs one, due only to the different scales shown. It is important 
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Figure 5: Same as in fig.^] but for T{Ao ^ bs + sb). 



to recall that there are some values of the widths that are not shown in the 
figures since they correspond to non allowed values for the squark masses, 
that is Mq < 150 GeV. In particular, in this first plot we can see that for 
the values of the parameters specified in the figure caption, values of about 
tan /3 > 40 are not allowed in order to keep the squarks masses in the allowed 
region. Keeping this in mind and the growing behaviour with tan/5, we can 
conclude that the searched parameter space region that maximizes the FC 
effect is localized at large tan (3 values, with the largest possible values being 
constrained to be below some maximum, which depends on the particular 
values of the other parameters. 

We next study the behaviour with rriA- We see again from plot (a) of each 
figure that it is clearly different, depending on the particular Higgs decay we 
are studying. The decay widths T{Ho —>■ bs + sb) and T{Ao — > bs + sb) clearly 
grow with niA due to the obvious phase space effects, that is, as uia increases, 
the corresponding Higgs mass increases too. In contrast, the T{ho — * bs + sb) 
width shows a less obvious behaviour with m^. As the lightest Higgs mass 
starts growing with and then it stabilizes, its decay width first grows with 
niA, then it reaches a maximum value and finally decreases with ttla, and we 
appreciate the setting of the decoupling behaviour for large tua values that 
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Figure 6: T{Ho tc + ci) in GeV as a function of (tan/5, (GeV)) 
(a), (fi (GeV), tan/5) (b), (^ (GeV), Mg {GeV)) (c) and (M^ (GeV^), A 
(GeV)) (d). The regions of the parameter space not plotted are the ones 
that give non allowed values for the squark masses (Mg < 150 GeV^). The 
values of the different MSSM parameters that have to be fixed in each plot 
have been chosen correspondingly to be: /i = —1500 GeV", Mo = 600 GeV^, 
Mg = 500 GeV, A = 800 GeV, tua = 250 GeV and tan/? = 10. 

will be better explained in section ^. 

The phenomenologically interesting values would be those that can 
allow the next planned colliders to detect and study all the three neutral 
Higgs bosons. For definiteness, in our numerical analysis whenever we have 
to fix niA, we have chosen = 250 GeV^, which for the interesting large 
tan P region, gives the three boson masses mh„, ttihoi and uia being accessible 
to the LHC. On the other hand, regarding particularly the lightest Higgs, as 
it is the one that has more possibilities of being detected in the next years, 
it would be also interesting to further study the region of lower uia values 
but we do not explore this here. 

Now we focus on the behaviour with ^. We can clearly see in figs.^ ^ and 
^ (b) and (c), that the widths are approximately symmetric under /i — > — 
and that for moderate < 600 GeV they all grow with this parameter. 
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Figure 7: Same as in fig.|] but for T{Ao — > tc + ct). 



We can appreciate this growing behaviour in all the flavour changing Higgs 
decays. The Ho and Ao decay widths grow with for all values, (figs.|^ 
and 1^ (b) and (c)); but the lightest Higgs decay first grows with from 
fi = 0, then, it reaches a maximum, decreases up to a minimum, and finally 
it continues growing (fig.^ (b) and (c)). We have found that this particular 
behaviour of the lightest Higgs boson is due to an accidental numerical can- 
cellation among the contributions from the different diagrams to the form 
factors, F^^{ho) and F^{ho), which takes place in this decay and not in the 
rest of decays under study. 

On the other hand, this behaviour with fi, as can be seen in figs. ^ and 
H (b)), depends also strongly on tan/5, because these two parameters appear 
together in the squark mixings, and the behaviour of the widths with them 
are correlated. For large tan jS they grow faster with fi, and for small tan (3 
the trilinear term starts to play a role too, so this growing smoothes down. 
Due to this correlation, one has to be careful in choosing independently both 
parameters as this could lead to non allowed values of the squark masses, 
Mg < 150 GeV. In particular, we can see again in plots (b), how for large 
values of the /i parameter, as for instance /i = 1500 GeV, only tan /3 values 
below about 40 are allowed. 
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Figure 8: Br{Ha ^ bs + sb), Ha = {ho, Ho, Ac) as a function of A for the 
selected MSSM parameters. Br{Ao — >■ bs + sb) is not plotted explicitly since 
in this range of the parameter space is undistinguisable from the Ho decay. 
Higher values of A are not allowed as they give Mg < 150 GeV. 



Following with Mg (figs.|^, ^ and ^ (c)), we find that the FC widths grow 
with Mg up to a certain value Mg^ and then they decrease very slowly (slow 
decoupling). Clearly, the value of Mg^ depends on the rest of parameters. For 
the explored MSSM parameter space region, we can see from these figures 
that it is in the range between about 200 GeV and 800 GeV. 

On the other hand, the behaviour with Mg is very clear too (figs.|^, ^ and 
^ (d)). All the decay widths decrease with this parameter, for large Mo, as 
it was expected. When Mo grows, the masses of the squarks inside the loops 
grow as well, thus the probability of producing such loop-induced processes 
decreases as they are suppressed by the squark masses in the propagators. 

To finish, we can see as well in figs, (d) that the partial widths in the d- 
sector are nearly independent of the trilinear parameters A for all the studied 
decays. 
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Figure 9: Br {Ha — > tc + ci), Ha = {Hg, Ao) as a function of A for the 
selected MSSM parameters. Higher values of A are not allowed as they give 
Mq < 150 GeV. 

In the following we will discuss the decay widths T(Ha — > tc + cf) (figs.|| 
and 0), where, as we have already explained, we will only study the behaviour 
of Ho and Ao decays, as the lightest Higgs cannot decay to tc or ct. 

Again, the first plot (a) of each figure represents the combined behaviour 
of the decay widths with tan/3 and m^. One of the fundamental differences 
among these decays of the w-sector and the ones studied before of the d- 
sector is that now there is a very mild dependence with tan (3 and the widths 
decrease softly with tan/?. Therefore, this time, the value that maximizes 
our effect is the smallest tan (3 value, but as we have said, the numerical effect 
of this parameter is not of crucial importance in these decays. 

Regarding the behaviour with uia, as can be seen in figs.H and |^ (a), it 
does not vary much with respect to the previous ones. The decay widths 
of the Ho and Ao again grow with uia due to the phase space as explained 
above, thus the same motivations to select a particular niA value will apply 
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here. 

The partial widths again grow with /i and are approximately symmet- 
ric under /i — —fi, but they show less symmetry than in the decays of 
the (i-quark sector studied earlier. The squark mixing term now goes with 
^t,c{^t.c ~ /icot (3) so the trilinear parameter A becomes more important for 
tan/? > 1. Looking at the plots | and ^ (b), (c) we can see that the max- 
imum values are reached for negative values of /i. We see as well that for 
large values of /i we can better appreciate the tan [3 dependence. 

Looking at the behaviour of the partial widths in the w-quark sector with 
the gluino mass, M^, in figs.^] and|^ (c), the widths again grow with Mg until 
they reach a maximum at a certain value Mg^, and then they start a slow 
decreasing (slow decoupling), as in the case of the ci-quark sector studied 
above. Again, the value of Mg^ also depends on the rest of parameters, and 
the mentioned behaviour is more evident for large values of the parameter /i. 

The behaviour with is again the same as before, decreasing with 
it as it was expected, but now, the trilinear parameter A acquires more 
importance than before, growing the flavour changing widths as we increase 
the value of A. Again here, it is worth mentioning that there are some 
values not being plotted that correspond to non allowed values for the squark 
masses. In particular, for small values of Mq, only small values of A are 
allowed. As a particular example, looking at figs.^ and |^ (d) we can see that 
for Mo = 400 GeV, values of the trilinear parameter A larger than about 
600 GeV would not be allowed. 

Finally, we study the behaviour of the corresponding FCHD branching 
ratios with A, and estimate their size for < A < 1 and for some selected 
points of the MSSM parameter space, which have been chosen belonging to 
the region where the PC effects we are looking for are maximized. This region 
can be extracted easily from our results in the previous figures. Thus, for 
the Ha —>■ bs, sb decays we select the following representative point: = 
250 GeV, tan/? = 35, /i = 1500 Gel/, Mg = 300 GeV, Mq = 600 GeV and 
A = 200 GeV; and for the Ha tc, ci decays we choose: = 250 GeV, 
tan/3 = 10, = -1500 GeV, M-g = 500 GeV, Mq = 600 GeV and A = 
800 GeV. Regarding the computation of the total Higgs boson widths, these 
have been evaluated with the HDECAY program | |37| |. 

For the Ha — >■ bs, sb decays, as can be seen in fig.|, the branching ratios 
grow with A, being exactly zero for A = 0, as expected, and can reach quite 
sizeable values, even for moderate A. For instance, for A ~ 0.6 the branching 
ratio for the ho is about 0.015, and for the Ao and Ho is about 0.25. Of 
course, larger A values would lead to larger rates, but these are not allowed 
here due to the Mg > 150 GeV restriction. 

For the Ha tc, ci decays, as can be seen in fig.|^, the branching ratios 
grow again with A, being exactly zero for A = 0. Again, only A values leading 
to Mq > 150 GeV are shown. The Ho branching ratio reaches its maximum 
value, of about 0.00035 for A ~ 0.8, while the pseudoscalar branching ratio 
reaches 0.0001 for this same value of A. 
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In summary, the FCHD branching ratios that we have found in this sec- 
tion are quite sizable, in fact, many orders of magnitude larger that the cor- 
responding SM rates (for instance, in the w-sector, and for ttihsm — 200 GeV 
we estimate B{Hsm ^ tc + ct) ~ 10~^^) and will produce an interesting 
amount of rare events at the planned next generation colliders. We do not 
study here this in more detail and postpone this interesting analysis for a 
future work. 



5 Non-Decoupling Behaviour of Heavy Squarks 
and Gluinos in Flavour Changing Higgs De- 
cays 

In the previous sections we have studied the FC effects of squarks and gluinos 
via SUSY-QCD radiative corrections in the decays of neutral Higgs particles 
within the MSSM. We have performed analytical and numerical calculations, 
and estimated the size of these effects as a function of the various MSSM 
parameters and A. In this section, we study the behaviour of these FCHD 
processes in the most pessimistic scenario where a heavy supersymmetric 
spectrum is considered, and we will show that these FC effects remain sizable 
even for squark and gluino masses as large as 0{lTeV). We will show that 
the reason for these corrections remaining so large is that at asymptotically 
large SUSY particle masses they manifest a non-decoupling behaviour. That 
is, the one loop SUSY-QCD radiative corrections to the neutral MSSM Higgs 
boson decays remain non-vanishing at asymptotically large squark and gluino 
masses, and it results in a finite and non-negligible contribution which in 
principle could provide an indirect signal of SUSY-QCD at low energies, 
even in this pessimistic scenario. 

The decoupling of all SUSY particles would imply that the prediction in 
the MSSM for all the observables involving non-SUSY particles in the exter- 
nal legs, such as the partial decay widths under study here, should tend, in the 
limit of large SUSY masses, to their corresponding values in the non-SUSY 
two Higgs doublet model, called 2HDMII in the literature . Formally, and 



following the lines of the Appelquist-Carazzone Theorem the decoupling 
would occur if the contributions of the SUSY particles to low-energy processes 
either fall as inverse powers of the SUSY mass parameters or can be absorbed 
into redefinitions of the couplings and parameters of the low-energy theory. 
In the present case, the decoupling of SUSY particles in the FCHD would im- 
ply, in particular, that the SUSY-QCD induced radiative corrections should 
tend to zero in the asymptotic large SUSY mass limit. This theorem has 
been proved to be valid for theories with an exact gauge symmetry, however, 
it does not apply to theories with spontaneously broken gauge symmetries, 
nor with chiral fermions, as it is clearly the case of the MSSM. Furthermore, 
in order to have decoupling, the dimensionless couplings should not grow 
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with the heavy masses. Otherwise, the mass suppression induced by the 
heavy-particle propagators can be compensated by the mass enhancement 
provided by the interaction vertices, with an overall non- vanishing effect, 
which is exactly what happens in some Higgs boson decays to fermions. This 
has been proved to happen in some ffavour preserving MSSM neutral Higgs 



boson decays in a series of previous works 1 40, 41, 42, 43, 44, 45] and also in 



production at hadronic colliders |^6|. Here we will show that the non 



decoupling effects also appear in the FCHD. These non-decoupling effects 
can have interesting phenomenological consequences as they imply that low- 
energy experiments can be sensitive to large mass scales, which cannot be 
kinematically accessed by direct searches. 

In order to show analitically the non-decoupling behaviour of squarks 
and gluinos in the FCHD, we perform a systematic expansion of the form 
factors involved, and hence of the partial widths, in inverse powers of the 
heavy SUSY masses and look for the first term in this expansion. To define 
our expansion parameter, we consider here the simplest hypothesis where 
all the soft-SUSY-breaking mass parameters and the fi parameter are all 
of the same order (collectively denoted by Ms) and much heavier than the 
electroweak scale. Mew, in such a way that the differences between these 
mass parameters are considered to be of order Mew- That is, 

M5 ~ Mo ~ Mg ^ ~ A > Mew (5.1) 

where again denotes the common soft breaking squark mass parameter. 
Mo = Mq = Mfj = Mf), and Mq, Mfy and Mjy are chosen here, for simplicity, 
to be the same for the two generations involved, and A is the common trilinear 
parameter. Notice that, with Ms ~ 0{lTeV), this choice will lead to a 
plausible situation where all the SUSY particles in the SUSY-QCD sector 
are much heavier than their SM partners.^] 

In the following, we perform the expansions of the SUSY-QCD contri- 
butions to the form factors, whose exact analytical results where presented 
in the previous section, in inverse powers of Ms and keep just the leading 
contribution of this expansion by considering that all the remaining involved 
mass scales mH^,mA,'m'h^,mz,mw and rUg are of order Mew- To this end, 
we use the results of the expansions of the one-loop functions and the rotation 
matrices that are given in Appendix B. 

For the Ha —>■ bs, sb decays, with Ha = ho, Ho, Ao, we find the following 
results for the leading terms in the expansions of the form factors, 

FYiho) = — tan /5 + cot aU^F A 

^ ^ ' 67r2mw/cos/3^ ^ ' M^ ^ ' 



^ Notice that the condition n ^ Ms is not necessary to get the SUSY-QCD particles 
heavy but is needed in order to get all the charginos and neutralinos heavy in the SUSY- 
electroweak sector. The condition of large trilinear couplings Aq ~ Ms is not necessary to 
get large SUSY masses. However it is a plausible choice from the theoretical perspective 
if one assumes that all the soft breaking parameters have a common origin. 
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T^hsz, X Gig nis sin a , ^ xU-^a-r^,^^ 

FR^ho) = -^-^ -tan/3 + cota^FA 

DVT 2mw cos p Mq 

F^'(Ho) = — tan/? -tana A 



FlHEo) = — tan/?-tanaU^F A 

^ ^ ' 67T 2mw cos ^ ' ^ ' 

Fl^{A,) = -z^-^tan/5(tan/5 + cot/5)^F(A) 



F'^{A^) = z^-!^tan/3(tan/3 + cot/3)^F(A) (5.2) 



Similarly, for the —>■ tc, ct decays, with Ha = Ho, Ag, we find the 
following results for the leading terms in the expansions of the form factors, 

Fl\Ho) = — — (cot/?-cota)^^F(A) 



Gtt 2mw sin /3 M^ 



Fj?(Ho) = — cot /? - cot a A 



= -^^-^cot/5(tan/5 + cot/5)^F(A) 
Fj,'{Ao) = ^^-^cot/3(tan/5 + cot/3)^F(A) (5.3) 



Some comments are in order. First, we can see from eqs.( |5.2| ) and ( |5.3|) 
that, taking all SUSY mass parameters arbitrarily large and of the same 
order, the SUSY-QCD contributions to the FCHD partial widths lead to 
a non-zero value. That is, they do not decouple in the large SUSY mass 
scenario. This can be seen clearly, for instance, in the simplest case of equal 
mass scales, /i = Mg = Mq, in which the dependence on the SUSY mass 
scale of the leading contributions disappear, and these provide a constant 
and non- vanishing value. 

In the above eqs.( |5.2D and ( [5.3[ ) there appear one new function that is 
defined as, 

F(A) = ^[(A + 1) ln(A + 1) + (A - 1) ln(l - A) - 2A], (5.4) 



where A is the parameter introduced in eqs.( p.l|) and ( p.2|) , and it is the 



only function carrying the information of the PC strength in our asymptotic 
results. Notice, that it gives a good approximation of the behaviour with 
A of the exact result for the FCHD branching ratios found in the previous 
section and shown in figs.^ and |^. For small A values, A <^ 1, it behaves as 
F(X) ~ — 2A/3 — A'^/5, and the first term, which is linear in A provides the 
corresponding result if the mass insertion approximation would have been 
used instead. 
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Figure 10: Non-decoupling behaviour ofT{Ha —>■ bs + sb) in GeV with Mq = 
H = A = Mg = Ms, for Ha = ho (top panel) and Ha = Ho, Ao (bottom panel) 
and for tan/3 = 35, A = 0.5, niA = rnn^ = 250 GeV and nih^ = 135 GeV. 
Exact one-loop results in solid lines for the ho and Ho, long-dashed lines 
for the Ao- The expansions given in eg. l \5. ^) (short- dashed) are plotted for 
comparison. 



The previous results of eqs .( 15.21 ) and (|5.3|) are valid for all and tan/5 



values and keep all the involved quark masses, m^, mb, rric 
from zero. 



and nis, different 



These expressions are of O 



■ Mew 
Ms 



)° in the large SUSY mass expansion 



and have corrections coming from the next to leading terms that are of order 
( ^if^ )" with n > 0. These latter are not shown explicitly since they vanish 
in the asymptotically large SUSY mass limit and therefore, they decouple. 

The results for the form factors in eq. ( p.2| ) are in agreement with those 
in the literature regarding Higgs-mediated FCNC effects, being radiatively 
induced from SUSY-QCD particles, within the context of B meson physics p, 
H, |l^ . These works use instead the effective Lagrangian approach, or, 
equivalently, the zero external momentum approximation, are focused on 
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Figure 11: Non- decoupling behaviour ofT{Ho ^ tc + ci) in GeV with Mq = 
fi = A = Mg = Ms and for tan (3 = 35, A = 0.5 and niHo = 250 GeV . Exact 
one-loop results (solid) and the expansions given in eg. ( \5. Sj ) (short-dashed 
lines) are plotted for comparison. T{Ao ^ tc -\- at) is not plotted explicitly 
since in this range of the parameter space is undistinguisable from the Ho 
decay. 




the large tan (3 limit, and use the mass insertion approximation, valid for 
small A <^ 1 values, to describe the flavour mixing of the internal squark 
propagators. Our results are more general in that they are for arbitrary 
tan/3 and A values, and they recover the previous results when the limits 
tan/? ^ 1 and A ^ 1 are considered in eg. (^-21). In addition, we include 
the predictions for the form factors which are proportional to and 
are usually neglected. In eg. ( p. 3D we also provide predictions for the form 
factors of the w-sector, Fj^ and Fj^, which are proportional to and rric, 
respectively, and are obviously not relevant for the 5-meson physics but can 
have interesting applications for future indirect SUSY searches at the top 
guark physics. 

It deserves special mention the convergence, in the » mz limit | i7[] , 
of our result for the ho form factors to the vanishing tree-level prediction of 
the SM rate. Actually, since cot a — tan/3 in this large niA limit, we can 
see this vanishing explicitly in Fl^j^{ho) of eg.(|5.2|). 

Finally, in order to show numerically this non-decoupling behaviour, we 
plot in flg]l^ the exact one-loop results for the Ha — >■ bs, sb partial widths 
and in flg.|lT] for the Ha tc, ct partial widths as functions of the large SUSY 
mass scale, Ms (solid lines for all except for Ao —>■ bs, sb that is plotted by 
long-dashed lines). For comparison, we also plot the approximate results 
from the leading terms of our expansions in egs.( |5.2| ) and ( [5.3| ) (short-dashed 
lines). The values of Ms have been taken here up to unrealistic very large 
values just to illustrate the convergence of the exact and approximate results. 
In all these plots, we choose Mq = Mg = fi = A = Ms and tan/? = 35, 
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A = 0.5. We also fix ruh^ = 135 GeV, mu, = 250 GeV and tua = 250 GeV. 
We can see in these figures that for large SUSY mass parameters, the exact 
partial widths tend to a non-zero value and can be quite large even for quite 
heavy SUSY particles. For example, T{ho — >■ bs + sb) (top panel in fig.p!0|) can 
be as large as 3 x 10~^ GeV for Ms — 10 TeV. The Ho/Ao — >■ 6s, sb cases are 
plotted together (bottom panel in fig.|10D as they are very similar, and can 
be as large as 7.6 x 10"^ GeV for Ms — 10 TeV . The same behaviour can 
be seen for T{Ho/Ao — > tc + cf), fig. 11, but these are smaller than in their bs 
decays, because of the phase space suppression and the unfavourable value 
of large tan/3. For example for Ms — 10Te\^ they are 1.9 x 10~^ GeV . 

In summary, the non- decoupling behaviour found in this section explains 
the large values for the FCHD partial widths obtained in the previous section 
for large values of the MSSM mass parameters. It is worth noticing that this 
non-decoupling property of the SUSY particles is associated to the fact that 
the mass suppression induced by the heavy sparticle propagator is being com- 
pensated because there is a Higgs-squark-squark coupling, eqs.( [A.lD , ( [A. 21 ), 
(|A.3| ) and (|A.4|) , with mass dimension that depends on /i and on A and there- 
fore grows with Ms in the large SUSY mass limit. On the other hand, from 
the numerical comparison between the exact and our approximate results of 
eqs.( ^.2D and (|5.3| ), we can conclude from fig.|lO]and fig.|ll| that these large Ms 
expansions, with just the leading term, are a good approximation for large 
enough SUSY mass parameters. These asymptotic results for the FCHD form 
factors can have interesting applications for future indirect SUSY searches 
at next generation colliders, which together with the previously found non- 
decoupling effects in flavour preserving Higgs decays |5D,5T,i^,^,^,^] and 



production 



could provide an indirect signal of SUSY, even in the 



most pessimistic scenario of a heavy SUSY spectrum at the TeV scale. 



6 Summary and conclusions 

Our main goal in this work has been to study the Flavour Changing Neutral 
Higgs Decays in the MSSM as a possible indirect search for supersymmetry. 
For that purpose we have looked for particular channels whose branching 
ratios were enhanced in supersymmetry with respect to their predictions 
in the SM and in the 2HDMII. If such processes are finally seen in any of 
the next generation colliders, they could provide interesting clues to physics 
beyond the SM and the 2HDMII. 

In order to discern between the MSSM and the non-SUSY 2HDMII, one 
has to search beyond tree level, where the contributions of SUSY particles, via 
radiative corrections, change the predictions of the observables. In particular, 
in this work we have focused on those possible differences at one loop level in 
the neutral MSSM Higgs boson FC decays into second and third generation 
quarks, namely: Ha bs, sb with Ha = ho, Ho, Ao and Ha tc, ct with 
Ha = Ho, Ao- These are induced mainly by loops of squarks and gluinos, as 
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they are the dominant ones since their size is governed by ag- 

In order to study these loop induced flavour changing effects, we assumed 
here the most general hypothesis of misalignment, where no extra symmetry 
has been included to simultaneously diagonalize the quark and squark sectors, 
leading to non-diagonal squark mass matrices. In this work we only keep 
the well motivated intergenerational mixing between cl and II for the up- 
type squarks, and sl and ht for the (iown-typc squarks, parametrized 'a 
la Shcr' by an unique parameter A which we take to be < A < 1. We 
performed the complete analytical calculation of the FC partial widths and 
studied numerically the size of these loop induced FCHD as a function of 
our simplified choice for the MSSM parameters, namely m^, tan/3, /x, Mg, 
Mo and A. In order to understand the behaviour of T[Ha — > fes + sh) and 
T{Ha ^tc + ct) in different regions of the MSSM parameter space, we have 
analyzed in full detail the dependence of such FC processes with all these 
MSSM parameters, being scanned by pairs, and for a fixed A value. After 
that, we have studied as well the behaviour with A, by varying it in the range 
< A < 1. 

From this numerical analysis we have learned that the hs decays grow 
with II and tan /3, being almost independent on the value of the parameter 
A, while the tc decays have a very mild dependence on tan/3, decreasing 
slowly with it, and the parameter A acquires more importance. Concerning 
the behaviour of the decay widths with we found that all of them decrease 
as Mo grows. Notice that we have been very careful when choosing different 
values for /x, tan /3, A and Mq since some choices do lead to too low values of 
the squark masses, being below our required value of 150 GeV. 

Looking at the behaviour with m^, we saw that due to the obvious phase 
space effect, the FC decay widths for Hg and A^ clearly grow with itla- In 
contrast, since the lightest Higgs mass has an upper limit, the FC decays 
involving ho start growing with niA but then decrease with it for large ruA 
values, manifesting the setting of the expected decoupling behaviour. On 
the other hand, the phenomenologically interesting value of would be 
one that can allow the next generation colliders to detect and study all the 
three neutral Higgs bosons. 

Finally, the behaviour of all the decay widths under study with Mg is 
clear, they first grow with this parameter, reach a maximum, and then they 
suffer a slow decoupling. 

With this complete analysis, we ended up by selecting some points of the 
MSSM parameter space that belong to the region where the these FC effects 
are maximized. Then we studied the behaviour of the corresponding branch- 
ing ratios with the parameter A in the previously selected region, obtaining 
quite sizable rates, whose implications at next generation colliders would be 
interesting to further study in the future. Concretely, for |A| < 0.6, and 
for the selected points in the MSSM parameter space specified in section 4, 
we have found the following approximate ratios Br{ho bs + sb) < 0.01, 
Br{Ho, Ao^bs + sb) < 0.2 and Br{Ho, Ao^tc + ci) < 0.00005. 
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Finally, we studied the behaviour of these FC decay processes in the 
most pessimistic scenario where a very heavy supersymmetric spectrum was 
considered and we showed that the size of these FC effects remain sizable even 
for squark and gluino masses as large as 0{lTeV). We have shown in full 
detail by an explicit analytical computation that this unexpected large size 
of the FC effects is due to the non-decouphng behaviour of the heavy squarks 
and gluinos in the loops. In fact, we have presented the asymptotical results 
corresponding to large SUSY masses, Ms >> Mew, for the FC form factors 
(and hence, the corresponding FC effective couplings, Habs and Hate), which 
can be of much interest for future estimates of rare processes production rates 
at the next generation colliders. Finally, we have studied the behaviour of 
T{ho ^ bs + sb) in the asymptotic limit of very large pseudoscalar mass, 
rriji >> rriz, and we have recovered, as expected, the vanishing SM tree-level 
result. 

In conclusion, the one loop induced flavour changing Higgs decays offer 
a promising scenario where to look for indirect signals of supersymmetry. 
These decays get quite sizable contributions in some regions of the MSSM 
parameter space, which remain, due to the non-decoupling behaviour of the 
squarks and gluinos, of considerable size even in the most pessimistic case 
of a very heavy SUSY spectrum. With such analysis in mind, and due 
to their interesting phenomenological implications, it would be crucial to 
further study the experimental possibilities of finding these effects in the 
next planned colliders. 
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Appendix A 



In this appendix we present the explicit values of the Higgs-squark-squark 
couplings in the squark mass eigenstate basis that appear in eqs.(|3.3|) and 
( |3.4| ). For the up-type squarks, these couplings follows: 



9HaU. 



aUfj 



gMz 

cos 6 
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l,a 



1 2 2 



My/ sin 13 
2Mw sin p 
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grric 
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with: 



l,a 
yu 

3, a 
I/" 



2Mvi/ sin P 



(sin(a + P), — cos(a + /?), 0) , 
(cos(a), sin(a;), 0) , 
{Aq cos(a;) + fi sin(a), Ag sin(Q;) 
{Ag cos(a) + /i sin(a), Ag sin(a) 



(A.l) 



— /icos(a), {Ag cos/3 + /isin/5)/z) , 

— /icos(a), —{Ag cos P + /isin/5)/i) 

(A.2) 



for Ha = {ho, Ho, Ao), respectively, and where Ag = At, Ac correspondingly. 
Similarly, for the down-type squarks, the couplings are as follows: 
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with: 



' l,a 



2, a 
d 

3, a 
d 

4, a 



- sin(a + P), cos(a + /?), 0) , 

- sin(Q;), cos(a), 0) , 

-Ag sin(a) — ji cos(a) , Ag cos(a) — /i sin(a) , {Ag sin /? + cos P)/i) , 
-Ag sin(a;) — ji cos(a;), Ag cos(a;) — /i sin(a), — {Ag sin/5 + yU cos P)/i) 

(A.4) 



for Ha = {ho. Ho, Ao), respectively, and where Ag = Af,, Ag correspondingly. 

In the previous formulas, i?^^ and with a, P = 1, 2, 3, 4, are the rota- 
tion matrices that diagonalize the squark squared mass matrices of eq.( p.l|) 
and (p.2|) respectively. 

Appendix B 

In this appendix we give the expressions required to compute the leading 
contribution to the FCHD partial widths in the large SUSY mass expansion 
defined in Sect. 5. For that purpose we first write the values of the squark 
masses and rotation matrices and then the formulae for the two- and three- 
point integrals. 
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The expressions for the squark masses and rotation matrices, in the hmit 
of large SUSY mass parameters and keeping just the leading contribution, 
are 

^M2(l + A), M| ^M2,M? ^M2, M| ^ M2(1 - A) (B.l) 
R[f ^ ^ 4f ^ 4? 7ffeF(^ - /itan/3) 

-0-14 — JTLii — rt23 — -"-22 -"-34 — -"-31 — -"-43 — -"-42 — ^ 

4J ^ ^ -4? ^ -4? ^ -7iiij(^ - /i tan/3) (B.2) 

with similar results for R^'^^ just replacing b ^ t, s ^ c and tan/5 cot /3. 

In this limit the expressions for the two- and three-point one-loop integrals 
involved are, 

CnK'<,mJ;M|,Mj,Mj) ^ ^/io(i?,., i?,J + O 

i?oK; Ml Ml).. A- log^ + g,{Rj + O (^) 
5iK; Ml,, M|) ^ -1a + + g,{Rj + O (^) , 

(B.3) 

where Rq^ = Mg/Mq^ and the explicit expressions for the functions fi and gt 
can be found in |^,^. For the simplest case, where Rq^ = Rq^ = 1, they 
are 1) = /io(l, 1) = /i3(l, 1) = 1 and g,il) = (72(1) = 0. 
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